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Abstract 

We provide an overview of Tsallis statistics, presented as a special case of superstatistics and applied to the mul¬ 
tiparticle processes described by the statistical cluster model. This model combines Boltzman statistics applied to 
hadronization of clusters and superstatistics induced by fluctuations of their Lorentz factors. It results in a Tsallis-like 
distribution of transverse momenta with some scale, the origin of which is discussed. We show that this distribution 
becomes a Tsallis distribution in a special case, namely when one combines fluctuations of the Lorentz factor and 
temperature inside the cluster, given by beta and gamma distributions, respectively. 
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1. Introduction 

Quasi-power like distributions, exemplified by the 
Tsallis distribution |[ll]. 


HE) = 


l-H-q) 


T\ 


( 1 ) 


are nowadays recognized as a legitimate tool for use in 
describing different distributions of this type (usually with 
E replaced by transverse momentum pj). This is the more 
so when it was realized that, albeit based on nonextensive 
statistical consideration^ it is equivalent to the so called 
Hagedorn distribution yj (known in the held of particle 
physicsfl 

n - 1 / £ 

HE)^-^[l.-] , ( 2 ) 


riT 


(and references thereiifl It is remarkable that such a dis¬ 
tribution also fits recent LHC data on transverse momen¬ 
tum distributions extending up to pr - 200 GeV, with 
cross sections dropping through ~ 14 orders of magni¬ 
tude, cf. for example ifnl] . There are numerous ways 
besides nonextensive statistics which lead to distributions 
o or 0, like, for example, fluctuations of the scale pa¬ 
rameter (especially the so called superstatistics), prefer¬ 
ential attachment or multiplicative noise - to name only a 
few (see 112, O] for recent reviews j§. 


2. Scale parameter fluctuations 

Let us concentrate on the case of superstatistics (under¬ 
stood as a superposition of two different statistics relevant 


where for n - both distribution coincide. The most re- 
cent examples of different phenomenological applications 
of one of these formulas can be found in [J, yl la [71 
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^One should mention that the first papers applying ^-statistics to scat¬ 
tering data in high energy physics dates back to @] and IToll . the second 
also pointing out the connection to Hagedorn’s theory. 

^Actually, the LHC results mentioned above were also interpreted 
in the spirit of a Tsallis-like distribution, emerging in this case, under 
some approximations, from the hai'd collisions of paifons described by 
QCD methods and followed by showering/hadronization processes into 
hadrons Ql. 
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to driven nonequilibrium systems with a stationary state 
and intensive parameter fluctuations) iQ. It is based on 
the observation that the Tsallis distribution ([TJ appears in 
a quite natural way as a result of T' scale fluctuations in 
the BG distribution, 

1 / E\ 

/(£) = -exp (3) 


when these fluctuations are describe d by a g amma distri¬ 
bution in the variable (l/T') |T^ 17, 18, 131. 
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nTT{n) \ T' / 


exp 


{- t ) 


q-V 


(4) 

where the parameter q defines the strength of the T' fluc¬ 
tuations, 

, , VariHV) 

Cj — \ - r —. (5) 

It should be mentioned at this point that depending on the 
statistical properties of the fluctuations, one obtains dif¬ 
ferent effective statistical mechanical descriptions. Tsallis 
statistics follow from the above gamma distribution of an 
intensive variable, but other classes of generalized statis¬ 
tics can be obtained as well and, for small variance of 
the fluctuations, they all behave in a universal way. From 
the thermal perspective, this corresponds to a situation 
in which the system considered, for example a decaying 
cluster, is not homogeneous, but has different tempera¬ 
tures in different parts, which are fluctuating around some 
mean temperature T (an essentially identical result is ob¬ 
tained when the temperature fluctuates from one cluster 
to another remaining constant in a given cluster). It must 
therefore be described by two parameters: a mean tem¬ 
perature T and the mean strength of the fluctuations de¬ 
fined by q. The gamma distribution (|4li itself emerges in a 
quite natural way, either from thermal considerations cor¬ 
responding to the purely multiplicative noise ifl^ . or in 
a stochastic approach without correlation between the as¬ 
sumed additive and multiplicative noises ra. 

We shall now discuss this situation in more detail. Note 
that the distribution g {T') in Eq. [U is in fact the product 


“^The first identification of local temperature fluctuations as the 
source of Tsallis distributions observed in high energy physics data was 
done in Gi, later on it was also used in analysis of data from cosmic 
rays events in IT^ . 


of two distributions: 


g{T')~gi{r)-g2{r), ( 6 ) 


where 

(r') = 

is a scale free power law, and 



g 2 {T') = exp 



(7) 

( 8 ) 


is an exponential, which for small values of T' cuts off the 
distribution g {T'), with the scale parameter T determin¬ 
ing how fast this cut-off is. 

Suppose now that we fluctuate the BG distribution (O 
using only the gi {T') part of the distribution g (T'). As a 
result we obtain the following scale free distribution: 

/ii(£) = ^ flfrgi(nexp(-^)cx£-^^'. (9) 

The scale appears when one cuts-off somehow the small 
values of T'. For example, a sharp cut-off of the small T', 
i.e., limiting T' to T' > T only, results in the following 
distribution: 


h2(E) = 


OC 



E-k+i 


F(y.-1)- 



where F(v, y) is an incomplete gamma function. The fac¬ 
tor F(/c - 1) - F(a' - l,EIT) now suppresses the power 
distribution (O for small values of E. However, the form 
of this suppression is not the same as in the Tsallis distri¬ 
bution. This can be seen by comparing the expansion 



r(^-i)-r(^-i,|) 


1 

- -H 
K 


^ y F(/ + ^) (EV 
^ r(/-i-/f-H i)F(/-h 1) \ t) 


( 11 ) 


with the corresponding expansion of the Tsallis factor: 



1 F(/ + k) IE V 

'i<^ F(i +K)r{i + i) r■ 


( 12 ) 


However, if we smooth out this suppression by replac¬ 
ing the previous sharp limitation of the integrand by some 
smooth suppression factor, provided, for example, by the 
exponential function g 2 (T') from Eq. (|8]l, we get, as re¬ 
sult, the Tsallis distribution defined in Eq. ([TJ. 
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3. Approximate description of cluster decay 


of two Bessel functions from Eq. (fTSl l. 


Recently, a new statistical cluster decay model of 
hadronization has been analyzed numerically, showing 
that even without resorting to approaches of the kind men¬ 
tioned above, the resulting distribution of transverse mo¬ 
menta follows rather closely a Tsallis distribution lEol] . 
albeit not identically so. In this work the hadronic clus¬ 
ters produced were supposed to decay purely thermally 
(i.e. by following the usual exponential Boltzmann-Gibbs 
(BG) statistics) but, at the same time, were supposed to 
move in the transverse direction with a fluctuating (trans¬ 
verse) Lorentz factor distributed according to the assumed 
power law. It turns out that the combination of both 
distributions follows (at least numerically) a quasi-power 
like distribution, closely resembling a Tsallis distribution. 
Note that, according to what was said before regarding 
superstatistics, production and decay of such clusters can 
be regarded as an example of superstatistics at work (and 
not necessarily resulting in a Tsallis distribution). 

In this work we shall investigate this phenomenon in 
more detail, aiming to obtain its analytical justification 
and a deeper understanding from the nonextensive statis¬ 
tical point of view. We start with statistical cluster decay, 
discussed recently in M- The distribution of transverse 
momenta proposed there ifl 

/((pr)oc J" dyrKo^yT^yoi^^rl - 

(13) 

where mr is the transverse mass and yj denotes the trans¬ 
verse Lorentz factor (Kq and Iq are modified Bessel func¬ 
tions of the, respectively, second and first kind). This fac¬ 
tor is then supposed to fluctuate according to some power¬ 
like distribution, 

g(yT)~yT‘‘, (14) 

i.e. according to the gi function from Eq. dTl) with T' 
replaced by yj- Numerical calculations in 1201] show that 
whereas for large pj one has 

f{pT)~p7-\ (15) 

there is a suppression for small values of pr- To describe 
its origin note that the pr distribution given by the product 


fipr) ~ ^ y)’ 

depends strongly on the transverse Lorentz factor of the 
decaying cluster, yj, 

exp(-p 7 -) for yj 


/(Pt) 


Pt 


for large yj. 


(17) 


Because: 

• Ko{x)lQ{y) ^ for.r,y » 


• yr- ^ 2 ^ for large yjS 

• nij — Pt for pr » ni, 

therefore we can rewrite Eq. (fOl ) in the following way: 


fiPr) 




Pry/ ‘ = 


r(x)-r 


(18) 


Note that because yr is always limited to values yr > I 
only, there is a natural cut-off in Eq. (fOT l. which makes 
it similar and comparable to Eq. (fTOl) discussed above. 
As a result, here the power-like distribution of pr is also 
suppressed in the region of small values of transverse mo¬ 
menta. The factor [r(/c) - T (k,ptIT)\ introduces the be¬ 
havior ~ pj for small pr and remains constant for large 
values of px- To summarize, such an approximation must 
therefore lead to a result which is very near to the numer¬ 
ical results presented in lEoll . and not coinciding with a 
Tsallis distribution. 

Because of the important role played by fluctuations of 
the Lorentz factor yj described by the distribution g (yj) 
in Eq. (fTST l. it is interesting to speculate on its possible 
origin emerging from the quark structure of the colliding 
nucleons. Let us denote hy x - prlph the fraction of 
transverse momentum of the parton, pr, with respect to 


^It is worth noting that cluster models have a very long tradition, a 
similar formula was already derived in I2II1 . 


®For yj > 5 the difference is less than 1%. We expand in series in 
z = l/yr: I/z - = 1/z - Vl - z^/z - z/2. 
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the momentum of hadron, ph, and let us assume that the 
density of the parton distribution is 


w(x) ^ Ax^(l - x). 


yr = 


Xl + X2 
2^jXlX2 


X\ - p cos < 


X2- p sin (, 


in which 

Ji :2 


= tant 


■*1 


: p^ and jT 


1 + tan (j) 

2 Vtan 0 


1 do¬ 
er djT 


= A 




1 + tan^ (j) 


‘ 1 + tan , tan 0 

1 “ P I = + P - 2 

Vl + tan2 (j) 1 + tan^ 

1 + tan c 


dXjT 


Denoting 


we can write 


y^jT 


2 Vtan 0 , 

1 + tan (p 

2 -\/tan0 


Jd<p(...)6(y) = Jdy(...)j^-^ (25) 




where po, are solutions of Eq. (l24l i for y = 0. Now, 


dy 


1 - tan (p 


It means therefore that Eq. (l2?t has now two terms with 
different limits of integrations over p in each of them. 


(19) 


In the center of mass system the Lorentz factor of the clus¬ 
ter formed by the collision of partons with fractions of 
momenta equal to, respectively, xi and X 2 , is given by 


pe(0,pi); Pi = 


1 


sin (poi 

(corresponding to ep e (nlA,nl2)) 
1 


-^1 -H tan^ 001 
tan 001 


( 20 ) 


p 6 (0,P2); P2 = —= Jl -Htan2 0o2- 

cos 002 '' 


(29) 


(30) 


Notice that if x\ - X 2 then jj - 1, it becomes greater 
than unity only for nonsymmetric values of xi, 2 - It is con¬ 
venient to change variables to 


( 21 ) 


■ ( 22 ) 


(corresponding to 0 6 ( 0 , 7 r/ 4 )). However, in both cases 
we cover the whole range of 2 Ci_ 2 . For the choice ( l29b it is 
done preserving all the time inequality X 2 > xi, whereas 
for the choice (l30l) it is done by keeping all the time 
X 2 < Xl. In both cases the final result must be the same. 
Therefore, choosing for example, Eq. we have that 


1 dcr 
(T djT 


= 2A^ 


1 


tan 002 

ftv I \ 1 H- tan^ 002 


d(l> 


a ^p-, 

Jo ' 


The Jacobian of this transformation is equal |7| = p. The 
distribution of the Lorentz factor jt is now given by 


<pQ2 


1 + tan 002 
1 -P , 

V Vl +tan2 0o2 


+ P 


dpp 


tan 002 

1 -I- tan^ 01 


l+2fl 


(31) 


For yj » 1 it can be numerically approximated by 


1 dcr 1 

cr djT -fj - 


^ Jt 


(32) 


(23) 


(24) 


where for a - I and b - 3 one has k - 5. These value 
seem to be a reasonable first guess in the case of nucleon- 
nucleon collisions (b - 2m- 1 , where m denotes the num¬ 
ber of spectator quarks, whereas the parameter a is given 
by the expected mean value (x) = 1/3). The exact form 
of Eq. (l32l i for a - I and b - 3 would be 


1 do¬ 
er djT 


with 


= 400- 


, ,_: (1 H- tan^ 0 ) (26) 

c/0 4 tan 0 ytan0 ' ' 

and there are two solutions, 0 o(i, 2 ), for which we have 

tan 00 ( 1 , 2 ) = 27 ^ - 1 + 2yT ^yj - 1, (27) 

where tan 0 oi > 1 and tan 0 o 2 < 1 . (28) 




: tan 002 


1 

35 


3r 

^ 42 




2 r^ 


• 2 rr V^ 




(33) 


(34) 
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4. Fluctuations of relativistic temperature 

We close these remarks with the following observation. 
In Eq. (fTSl l one encounters the scale factor Tyr, not yr 




































alone. Some time ago such a quantity was proposed as 
relativistic temperature 0221] . 


r* = jtT, 


(35) 


(albeit in this case it would rather be a relativistic trans¬ 
verse temperature, whatever that means). From the dis¬ 
cussion following Eq. (fTOl i it follows that to get a proper 
Tsallis distribution when fluctuating the scale parameter, 
one has to use the full gamma function, Eq. dUi. There¬ 
fore, fluctuating Tyj should bring the distribution of the 
transverse momenta to the desired final Tsallis form. This, 
however, would mean fluctuation of the relativistic tem¬ 
perature, T*, but so far its proper form is still under dis¬ 
cussion 1231]. Nevertheless, keeping this reservation in 
mind let us elaborate on such a possibility in more detail. 

Assuming a beta distribution for fluctuations of 1 jyj. 


1 


1 


-f- 

jt} r(A--I-l)r(a h- 1) (yj- 


1 


1 

Tr 


and a gamma distribution for fluctuations of l/T, 

\t) too \ t } T r 

one gets the joint distribution 




(36) 


(37) 


(38) 


This means that one could therefore consider the existence 
of fluctuations of the relativistic temperature l/T*, which 
are again given by a gamma distribution with parameter 
K defining the size of these fluctuations. Relative fluctua¬ 
tions, defined as 


ojiz) = 


Var(llz) 
<l/z)2 ’ 


(41) 


satisfy the relation 


w(n 


tj(T) + oj (yr) + 2oj{T)(jj (jj) 
1 H- 3w(7’) 


(42) 


which connects fluctuations of T and yj 0 . The corre¬ 
sponding nonextensivity parameter is 


1 = -=a;(r). 

n 


(43) 


Note that the parameter n, which defines fluctuations of 
temperature in the decaying cluster, cf. Eq. (|4|, is pro¬ 
portional to heat capacity of the system under constant 
volume, n oc Cy. In the case considered here it is much 
greater than k, n » k. On the other hand, as shown in 


I 2 LI 25 II . the observed power index in the corresponding 


Tsallis distribution is smaller than n. The above consid¬ 
eration provides a simple and natural explanation of this 
fact, the power index is defined by both n and k (i.e. by 
fluctuations of the relativistic temperature). 


and can define the function 

^ (nTof 

T*^ yj) r(n)r(/f-H l)r(a) ' 



In the case when the parameters of the two components 
of this joint distribution are related in a certain way, for 
example \f a - n - \ - k, fluctuations of the relativistic 
temperature T* are again given by a gamma distribution 
(but this time with a changed shape parameter): 



5. Concluding remarks 

To summarize, the statistical cluster decay mechanism 
proposed in 020(1 and discussed above, is yet another 
example of supertatistics (which is not necessarily con¬ 
nected with a Tsallis distribution). Eluctuations of the 
Lorentz factor given by the distribution g {jt) fully spec¬ 
ify the slope parameter of the transverse momentum dis¬ 
tribution in the region of large values of pj. If they are 
given by a scale free power-like distribution, as in Eq. 
Gl, the resultant distribution of pj is also a scale-free, 
power-like one. Its behavior for small values of pj is dic¬ 
tated by the fact that the Lorentz factor is defined only 


^This is another example of the generalized thermodynamic fluctua¬ 
tion relations discussed previously in l24h . 
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for jT > I, therefore there is always a natural cut-off in 
g {jt), eliminating yj- < 1. As a result one gets a distribu¬ 
tion which is not a Tsallis distribution, remaining, how¬ 
ever, quite close to it numerically. 

One can therefore invent a more general distribution, 
which is a product of a beta distribution in l/yj- and a 
gamma distribution for the parameter IjT. If the parame¬ 
ters dehning these distribution are related in some specihc 
way, the resultant distribution is again a gamma distribu¬ 
tion with a modihed shape. In this case one gets the true 
Tsallis distribution for the spectrum of the transverse mo¬ 
menta with the slope parameter determined by n* from 
Eq. (l43l l. Note that in this case it can be argued that it 
is determined by fluctuations of the relativistic tempera¬ 
ture T* = yjT (modulo reservations concerning its proper 
form as mentioned before). Combining this result with 
the connection between fluctuations of temperature and 
specihc heat, (x>{T) ^ 1/Cv, one realizes that huctuations 
(jj {yj) provide a natural source of Huctuations resulting 
in the observed values of the slope parameter n*, which is 
smaller than n ^ Cy. In such a way the puzzle discussed 
in 13, 12, 23] could Hnd its natural solution. 


We close by concluding that, as shown here, the real¬ 
istic superstatistics for high energy physics could have 
deviations from the standard Gamma distribution and 
hence from Tsallis statistics, although in an approximate 
sense Tsallis statistics is still recovered. This provides yet 
another support for the idea expressed in ll^ that general 
superstatistics, based on more general distributions than 
Gamma distribution, might be relevant in high energy 
physics. 
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